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^ : ABSTRACT 

O ! The structure and magnitude of the electric field created by a rotating accretion disk with a 

O [ poloidal magnetic field is found for the case of a vacuum approximation along the axis. The 

accretion disk is modeled as a torus filled with plasma and the frozen-in magnetic field. The 
dimensions and location of the maximum electric field are found, as well as the energy of the 
^ , accelerated particles. The gravitational field is assumed to be weak. 

> 

)0 : 1 INTRODUCTION 

ON 

' Recently, there has been wide discussion of various mechanisms for accelerating particles around 
supermassive black holes (SMBHs) in the nuclei of galaxies and stellar-mass black holes in our 
Galaxy, in connection with studies of synchrotron radiation and inverse Compton scattering 
in the well collimated jets observed from radio to gamma-ray wavelengths. Extremely high- 
angular-resolution observations obtained via radio interferometry show that these jets become 
very narrow (comparable to the gravitational radius) with approach to the black hole. Explana- 
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c/3 . tions of particle acceleration near relativistic objects (black holes and neutron stars) are usually 
^ ' based on two types of mechanisms: acceleration by electric fields and magnetohydrodynamical 
acceleration (the Blandford-Znajek mechanism [1]). Acceleration by an electric field, and the 
^ . very existence of the electric field, are inseparably linked with the low density of plasma in this 
volume. Conditions justifying a vacuum approximation are probably realized in the magneto- 
spheres of pulsars and in some types of SMBHs [2, 3]. In this case, it is possible to accelerate 
particles to extremely high energies [4]. The limiting charge densities for which the vacuum 
approximation remains valid are determined by the formula [5] 

(l/p).(i//lo*)■lo-2-^ (1) 



27rce 

Here, fl is the angular velocity of rotation, H the characteristic magnetic-field strength, P the 
rotational period, c the velocity of light, and e the electron charge. It is evident from this 
expression that we should have for typical quasars rig < 10~^cm~^. Note that, in intergalactic 
space, rif. ^ 10~^cm~^, and in the Galaxy, ^ Icm^^. The presence of a black hole in the 
center of the Galaxy also leads to a decrease in rig near the center. In addition, the magnetic 
fields near SMBHs can reach values of the order of lO^G [4]. In any case, the question of the 
applicability of the vacuum approximation is rather complex, and must be solved taking into 
account the physics of black holes. In this paper, we will assume that the conditions for the 
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vacuum approximation are satisfied; this will enable us to investigate the structure of the electric 
field excited by a rotating accretion disk with a poloidal magnetic field. The formulation of this 
problem is analogous to that considered by Deutsch [6] , who found the structure of the electric 
field created by a dipolar magnetic field frozen in a rotating star. If a conductor rotates together 
with a frozen-in magnetic field, then, in a rotating coordinate system in which the conductor 
is at rest, there must be no electric field inside the conductor. Therefore, in an inertial system, 
an electric field is induced inside the conductor due to the presence of the magnetic field, and 
this electric field gives rise to a surface charge (in the special case of a magnetic dipole with a 
quadrupolar distribution). This surface charge is the source of the external electric field. We 
will consider the analogous problem for an accretion disk. 



2 CONSTRUCTION OF THE MODEL 

Let us consider a stationary system consisting of a rotating accretion disk in the form of a 
regular torus filled with plasma and surrounded by a poloidal magnetic field. We neglect 
motion of the torus toward the center. The geometry of the torus is determined by the two 
radii a and R (Fig. 1). Due to the high conductivity of the plasma, the magnetic field is frozen 
inside the torus, and currents flow only in the toroidal direction; there is no matter outside the 
torus. In a coordinate system comoving with the plasma, the electric fleld vanishes due to the 
relation 

fa = ^'^aO, (2) 

where j'^ = are the poloidal components of the current density, cr' is the plasma conductivity, 
and are the covariant components of the electric field tensor in a system comoving with the 
plasma. We introduce the following coordinates (Fig. 1): x is the distance from an arbitrary 
point to the center inside the torus in the same meridional plane, ip is the angle between the 
direction toward the center of the system and the direction to a given point from the center 
inside the torus in the same meridional plane, and (/? is the position angle defining this plane. 
Thus, the differential coordinates dx, dijj, and d(p form a righthanded orthogonal vector triad. 
We can write the square of a hnear element in Minkowski space for the differentials of these 
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coordinates ^: 

f ds"^ — dt^ — dx^ — x"^ dip"^ — {R — X cosip)"^ dip"^ , , 
1 = x\R — xcostpl- ^ ' 

Let the plasma in the torus rotate with angular velocity Vt with respect to a distant observer. 
Then, the coordinate transformation is given by ^ 

dx'^dx'^[5l + n5l5l] (4) 

We now introduce the covariant four- vector potential of the electromagnetic (EM) field Ai. Due 
to the axial symmetry of the system, only Aq, the potential of the electric field, and A^, the 
potential of the magnetic field, differ from zero. Therefore, the EM-field tensor Fjj has only 
poloidal components: 

FaO = daA^ , Fa^p = OaA^. (5) 

These components transform in accordance with (4) [7, §83]: 

F^Q^F^o + ^F^^, F^^^Fa^, A'o^Ao + nA^, A'^^A^. (6) 
Since F^q = in the plasma [see (2)], we have for x < a: 

FaO = —^Faip . (7) 

It follows from (5) and (7) that, when x < a, Aq = const — QA^p. It follows from the third 
equation of (6) that this constant is A^ inside the torus. 

The continuous boundary conditions for the tangential electric and normal magnetic com- 
ponents of the EM field act at the interface between the plasma and vacuum (at the surface of 
the torus). These components should vanish at the equator, due to the axial symmetry of the 
system and the mirror (anti-)symmetry of the components of the EM field. Hence, the normal 
component of the magnetic field can be expanded in a Fourier sine series: 

F^^ = nRn{x) sin(V'n). (8) 

Here and below, the summation over n is assumed, where n runs through all numbers of the 
natural series. It follows from (7) and (8) that the boundary condition for the tangential electric 
field is 

Ffo{a,^) = -flnRn{a) sm{i/jn) . (9) 



3 POTENTIAL AND STRUCTURE OF THE ELECTRIC FIELD 
NEAR THE TORUS 

Near the torus [see formulas (30), (32), and (37) in the Appendix], the main contribution to 
the potential Aq is made by the first term of the first harmonic of the Fourier series. Far from 
the torus, the potential dies away. The kinetic energy of a charged particle accelerated by the 
system is determined primarily by this part of the potential. Accordingly, we obtain the main 
approximation for the difference in the potentials at the torus surface between the angles 
and Ip — Ip: 

AAo ^n-R-Ho-a- [ln(4/6) - 1] • [1 - cos V-J/tt . (10) 

^We take the velocity of light to be unity: c — 1. 

^Where not indicated otherwise, Xi = t, x, tp, (p, nd Xa = x, ip. 
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Figure 2: Appearance of the electromagnetic field lines in the system. The solid curves show the 
magnetic field lines and the dashed curves the electric field lines, z and p are expressed dashed curves 
the electric field lines, z and p are expressed in fractions of R. 



This same expression can also be obtain in another way [8, §63]. The corresponding invariant 
result has the form 



e = WAA, = J [U' Fji] dx' = J U"'F^^ dip 
c 



Here, W denotes the four- velocity of the observer ^ at the measurement point, and the contour 
for the integral is chosen for convenience to be on the inside surface of the torus in a reference 
system comoving with the torus. After substituting into this expression If^ ~ Q, and the 
expression for F^^ = F^^p we obtain (10). 

The electric field inside the torus is zero, since the torus is conducting. In an incrtial 
reference system, this is accomplished by the compensation of two fields: that induced by the 
rotation of the magnetic field and the field of the surface charge on the torus. The surface 
density of the electric charge on the torus is Pe = ^-^a;o/(4vr). F^o is the normal component 
of the electric field on the torus surface [see (7)]. At large distances from the torus surface 
{x » a), these charges represent a superposition of dipoles (Fig. 2). However, at distances 
X » R, the field from all these dipoles has a quadrupole character. It is not difficult to obtain 
an expression for the dipole moment per unit angle (p: 



<i=^^^l2-ln(4/(,)l. 



(12) 



The electric field and its potential (p can be obtained by integrating all the dipoles over the 
angle (p [7, "40]. As a result, we obtain for the potential and components of the electric field 
in cylindrical coordinates the quadrature expressions 



0(p, Z) 



2d 



p cos if 



[1 + p^ + P — 2pz cos 93)2/2 



dip . 



(13) 



'- 



cos ip{z^ — 2^"^ — 2) + 3p + p cos^ (p 
{1+ + P — 2pzcos(py/'^ 



dip . 



(14) 



^Who is at rest relative to the distant stars. 
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Figure 3: Magnitude of the poloidal component of the electric field E^\{p/R, z/R) in units of its value 
at the saddle point. Z and p are expressed in fractions of R. 



Here, p = p/R and z = z/R are dimensionless cylindrical coordinates. The last three 
formulas are expressed in terms of derivatives of the full elliptical integrals. The result of this 
integration is shown in Fig. 2. 

Acceleration by the electric field is associated only with the component parallel to the 
magnetic field. This acceleration is especially efficient on the fl axis, where the electric and 
magnetic lines of force are parallel. We can readily see from (15) that the maximum electric 
field on the Q axis is reached at the point z = R/2. However, this point is not a local maximum 
of the modulus of the longitudinal electric-field component £"11 = (EH)/i7. Figure 3 shows the 
surface of the E\\ force field, where we can see that the point z = R/2 on the axis is a saddle 
point of the distribution of £"11. With approach to the torus, the conductivity of the medium 
should increase and the field should become force-free: E\\ 0. ^ 

According to (13), the potential at the saddle point is approximately 28% lower than the 
potential at the center of the system, and is roughly an order of magnitude lower than the 
maximum potential on the torus surface. The extent of the region of acceleration along the 
z axis (at the level of 0.5 of the value of E\\ at the saddle point) is determined by the points 
Zi ~ 0.25, Z2 ~ 1.25 in fractions of R. Assuming that a mass M with its corresponding 

gravitational radius rg is at the center of the system, we can use (13) to estimate the energy 
to which particles with the elementary charge initially at the saddle point on the axis can be 
accelerated: 

We can see that the factors in parantheses can be of the order of unity in quasars, so that the 
kinetic energy of particles accelerated by such a system can reach 10^° eV. 

*Note that, in a Blandford-Znajek model, the field is forcefree everywhere by definition, so that efficient acceleration 
is not possible. 



3z{l — pcosif) 



(1 + + — 2pz cos 1^)^/2 



d(p . 



(15) 
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4 CONCLUSIONS 



We can draw the following conclusions from our results. 

1) The magnetic field at large distances approaches a dipolar field, while the electric field 
corresponds to a quadrupolar distribution for the charge induced on the torus surface. It follows 
from Fig. 2 that, in the model considered here, in contrast to a Blandford-Znajek model, we 
find a tendency for the electric field lines to become more concentrated near the Q axis, which 
can explain the observed focusing (collimation) of the accelerated rclativistic particles. 

2) The dimensions and locations of the regions of cosmic-ray acceleration we have found can 
be used to compare our results with observational data. 

3) Thanks to its covariance, our method for the computation of the electromagnetic field 
in a system with toroidal symmetry can be generalized to the case of a gravitationally curved 
space-time. 

4) The mechanism considered here yields accelerated- particle energies with the same order 
of magnitude as the Blandford-Znajek mechanism (see, for example, [1, 9, 10]). 
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6 APPENDICES 

Finding the External Solution 

Let us write Maxwell's equations for arbitrary curvilinear coordinates in the axially sym- 
metric and stationary case [7, §90] : 

e'^^'^dpF^^ = ; dp {^rgg'^^g'^'^F^^) = An^r ; 
e-P^df,F^o^O; dfs {V^g^^g'^F^o) ^ f . ^ ' 

Here, e"^^^ is a Levi-Civita symbol, and g^^ are defined by expression (3), and f is the 

current 4-vector, which is identically equal to zero when x > a. We obtain for the magnetic 
field from (8) and the first of equations (17): 

Fx,p = -da:Rn{x) cos{ipn) , = -Rnix) cos{ipn). (18) 

Using (9), the external solution for the electric field can be expanded in a Fourier series in the 
variable ip. Then, in accordance with the third equation of (17), we obtain the solution for the 
electric field (for x > a) 

F^o = nZn{x) sin(^n), F^o = -dr,Zn{x) cos(^n), = -Zn{x) cos(V'n). (19) 

It follows from (8) and (19) that the external electric field is generated by the normal component 
of the magnetic field at the torus boundary. Everywhere where the expression {R — xcosip) 
is positive, such as in the range a < x < R, we can remove the modulus signs in the second 
expression of (3). With this in mind, introducing the dimensionless variable y — x/{2R) and 
denoting a derivative with respect to this variable with a prime, we substitute (19) into the 
fourth Maxwell equation (17) and obtain 

y [y^Z'l + 2yZ'^ - n{n 1)Z„] cos{(n - l)^'}- 

- [y^Z':^ + yZ; - n^Zj co^{^m]+ (20) 
+y [y'^Z'l + 2yZ'^ - n(n + l)Z^] cos{(n + 1)^)} = . 



6 



00 



I- - 



Figure 4: Contour passing around the system, shown by the dashed hne. 



The boundary conditions for Zn{y) when y — b = a/ {2R) and n > follow from (9): 

Zn{b) = -QRnib) , (n > 0) . (21) 
The necessary condition for the total electric charge in the torus to be zero can be written 

Q = y^i^xO^°V V^'^T^dip = {x : a<x<R). (22) 

This expression flows from integration of the fourth equation of (17). We thus obtain 

Z'o = yz[ (23) 

Setting the expressions with the same harmonics in (20) equal to zero, we obtain a recurrent 
system of equations for Zn{y)- Equation (23) is contained in the Maxwell equation (20) with 
harmonic n = 0. The condition that the total charge of the system be equal to zero can also 
be written for x > R. For this, we must take a contour passing around the system, so that 
{R — xcosip) does not change sign; i.e., the contour should be located to the right of the 
halfplane from the fl axis. We choose this contour as shown in Fig. 4. Thus, we pass by the 
torus along a semicircle with radius x > a on the right-hand side, and we allow the contour to 
approach infinity at the angles ijj — 7r/2 and ijj — — 7r/2 along half-lines parallel and antiparallel 
to the Q axis. The two (upper and lower) "caps" of the finite area vri?^ remain unclosed at 
infinity, but since the field there asymptotically approaches zero, the integrals on these "caps" 
can be neglected. We can write the surface integral over this contour in accordance with the 
fourth of equations (17) taking into account the mirror symmetry for the electric field in the 
equatorial plane: 



7i LXJ 



V'=7r/2 



0. 



(24) 



tt/2 X 

After straightforward but cumbersome manipulations, we obtain 



yz'n 



cos(^) sin(^) 
•n? — 1 n 



'TTn\ r Zn , 
+ sm|-)/-d, = 0. 



(25) 
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Let us solve this problem with accuracy to within the first three terms in the Fourier expansion 
{Zq, Zi, and Z2). In this case, taking into account (23), the last expression becomes 



7 Z 2 

l-2y')yZ[ + J -J-dy^-y'Z',. (26) 



Differentiating, we obtain 



(1 - 2y^)y^Z'; + (1 - 6y')yZ[ - Z, ^ \y{y^Z'i + 2|/Z^) . (27) 

We now write the Maxwell equation (20) corresponding to the harmonic n—1: 

(1 - 2y^)y''Zl + (1 - Qy^)yZ[ - Z^ - y{y'Z'^ + 2yZ!, - 2Z-,) = . (28) 

An equation for Z2 follows from these last two expressions: 

Z/X + 2yZ^- 6^2 = 0. (29) 

The solution vanishing at infinity has the form 

Z2{y)^C2ly\ (30) 

The constant C2 is found from the boundary conditions (21). We can find ^i(y) using any of 
equations (27) or (28): 

(1 - 2y2)y2Zf + (1 - ^y^)yZ[ -Z,^ ^C^/y' . (31) 

A partial solution of the inhomogcneous equation (31) has the form Z\ = 4C2/ (3y^). However, 

it does not satisfy (26), of which (27) is a consequence, so that we must search for a solution of 

00 

Ziiy) from (26). The substitution / {Z"/y) dy = f^o — f{y) brings the homogeneous equation 

y 

(26) into the form (1 — 2y'^)y{yf')' — / + /oo = 0. Hence, in the limit y — > 0, we obtain 
Zl - C,y^\ 

In the limit y — > 00, the substitution y — exp(^) reduces (26) to 

r(O+e"'^//2 = 0, r(0 + e-^«//2 = e-2«/oo/2, / = / + /. 

In the main approximation in l/y, the solution has the form f{y) = exp[— l/(8t/^)] + l/{8y^). 
Hence, the asymptotic for |/ — 00 has the form Z^ — > —Ci/{32y'^). 

Zo{y) can be found using (23) after Zi{y) is found. We present the asymptotics of these 
solutions for small and large y: 

^0 ^ -^1 + ^C2/{Sy) , Zi — ^ CJy + iC^/iSy^) , 
^o^-(Ci + 8C2)/(2%3), Z,^-{C, + 8C2)/{32y^). ^^^^ 

2/— >oo J/— >oo 

Here, Ci is the coefficient for the solution of the homogeneous equation (26); like C2, it is 
determined by the specific configuration of the magnetic field from the boundary conditions. 



Calculation of Coef ficients for the Special Case 
of a Current Ring in the Torus 
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Let us now find the magnetic field in a simple case. Let the current in the torus be in the 
form of a ring along the torus axis with a delta-function distribution [7, §90]: 

f = hm J^S{x - xo) S{ij - ^o)/(27r^) , = . (33) 

X0—^0 

To find the potential A^p corresponding to this current, we introduce the physical components 
of vectors in accordance with the definition 



Hphys = H = {-f^^v^ba/?!}- 



This is necessary because the Biot-Savart law in its usual form [7, §43] is written in physical 
coordinates: 



Here, e° is a unit vector in the direction of the angle (p, (j'T-e^) = j^\l\g^^\ cos(/9, and |r — = 

2R{R — X cos (f) + x^ is the square of the distance from the segment of current with coordinate 
(f to the observation point. 

Introducing the notation (p — 7r-|-20, — (1 — 2y cos'0)/(l — 2y cos^'+y^), we can transform 
(34) to the form 

f 2 sin — 1 

l-2ycosV'-K- / , ^ =d(f) 



"2 sin^ ( 



or 



= 2{rR/-K) ■ ^l-2ycos%l) ■ {k{k){1 - /2) - E{k)] /k . (35) 

Here, K{k) and E{k) are full elliptical integrals. 

We can find the asymptotic of (35) as y — > 0, or equivalently as — > 1. Further, using (18), 
we can obtain for the Fourier coefficients of the magnetic field: 



lim : - ^ [2 - ln(^)] ^R,^-^[y{l- ln(p}] ^R,^-^[y-{l- ln(^)/4}" . 

(36) 

It is clear from the boundary conditions (21) that, according to (36), \Z2{h) / Zi{h)\ b 
as 6 — > 0, so that we can neglect the remaining terms of the Fourier series in the case of this 
magnetic- field configuration. We introduce the notation Hq = \dxA'^\{x=R,i,=o} — J'^/{2R). 
for the magnetic-field strength at the center of the system. We find the coefficients in this 
approximation from (21), (32), and (36): 

Ci ^ 2nH^R%'^ [1 - ln(4/6)] /tt , 2nH^R%^ [1 - ln(4/6) /4] /tt . (37) 
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